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ABSTRACT 

Asymptotic Giant Branch (AGB) stars generally have spherically symmetric 
envelopes, whereas most post- AGB stars and Planetary Nebulae (PNe) show ax- 
isymmetric circumstellar envelopes. While various mechanisms for axisymmetric 
circumstellar structures may explain the shapes of PNe, they do not address how 
the shape of the circumstellar shell evolves. Here we address the temporal changes 
in the axisymmetry of AGB star envelopes, and in particular the development of 
the torus required in the Generalized Interacting Stellar Winds (GISW) model. 
Assuming (1) an AGB star rotates with sufficient angular speed at the start of 
the AGB phase; and (2) that the rotational angular momentum of the AGB star 
is conserved, we demonstrate that some very important observational features of 
AGB star axisymmetry evolution can be reproduced. We find that, compared to 
the star's increasing luminosity and decreasing effective temperature, the decreas- 
ing mass of the star primarily affects the axisymmetry of the envelope. When a 
representative mass loss history is adopted, where most of the mass is lost near 
the end of the AGB phase, the envelope's axisymmetry increases over time, with 
the strongest increase occurring near the end of the AGB phase. This may nat- 
urally explain why most AGB stars have spherically symmetric envelopes, while 
axisymmetry seems common-place in the post-AGB/PNe phase. The degree of 
axisymmetry at the end of the AGB phase is found to increase with increasing 
main sequence mass, and the onset of axisymmetry occurs only after the onset 
of the superwind (SW) phase, in good agreement with the observations. 



Subject headings: stars: AGB and post-AGB — stars: circumstellar matter - 
stars — evolution — stars: late-type — stars: mass loss — stars: winds, outflows 
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1. Introduction 

Asymptotic Giant Branch (AGB) stars are the highly evolved descendants of 1 ~M 
~8M main-sequence stars (Iben & Renzini 1983), which lose mass at a high rate (10~ 7 to 
10~ 4 M Q /yr; see e.g. Habing 1996; Habing and Olofsson 2004). The outflowing matter creates 
a dusty molecular circumstellar envelope which may completely obscure the central star. As 
the dust is driven outwards by radiation pressure of the central star, it drags the molecules 
in the envelope along with it. It is believed that mass- loss rate increases significantly only 
towards the end of the AGB phase (e.g. Vassiliadis and Wood 1993; Villaver et al. 2002a). 
This increase in mass-loss rate is necessary to explain the densities seen in typical Planetary 
Nebulae (PNe; Renzini 1981). The sudden and rapid increase in mass loss has been dubbed 
the superwind (SW; Renzini 1981). Since the invocation of the superwind, many observations 
of AGB stars have supported this hypothesis (e.g. Knapp & Morris 1985; Wood et al. 1992). 
During the SW phase the mass-loss rate, M, exceeds Msw ~ M Q /yr (e.g. Renzini 1981; 
Wood et al. 1992; Groenewegen and Marigo 2004). After the SW phase AGB stars enter 
the post- AGB phase, or pre-planetary nebula (PPN) phase, via which they may eventually 
evolve into PNe. During the post-AGB/PN phase heavy mass loss stops, and the envelope 
moves away from the star. Meanwhile, the central star heats up and ionizes the circumstellar 
gas. The remains of the envelope eventually disperse into the Interstellar Medium (ISM), 
while the central star cools off as a white dwarf (WD). 

AGB star evolution is intimately linked to its dust production. The circumstellar dust 
drives mass loss and shapes the circumstellar environments. While AGB stars generally 
have spherically symmetric mass loss (e.g. Olofsson 2004), most post- AGB stars and PNe 
show axisymmetric circumstellar envelopes (Meixner et al. 1999; Waelkens and Waters 2004; 
Sahai 2004; Soker and Subag 2005). The axisymmetry is seen as a bipolar or hourglass- 
shaped reflection nebula, where the two lobes (bubbles, jets) of the nebula are separated by 
an equatorial waist, or torus. The onset of the axisymmetry is believed to occur very near 
the end of the AGB phase (e.g. Olofsson 2004; Waelkens and Waters 2004). While some 
PNe are indeed spherical (Soker and Subag 2005), these appear to be old nebulae, and it 
is increasingly evident that the majority of PNe are aspherical, and that most are bipolar 
(Soker and Subag 2005). 

A suggested mechanism for shaping PNe is the Generalized Interacting Stellar Winds 
(GISW) model (e.g. Kwok 1982; Kahn and West 1985). In this model a circumstellar torus 
is present in the equatorial plane of the AGB star. The exact origin of the torus is not 
specified, but is agreed to be somehow created by the slowly expanding AGB wind. Later, 
during the post- AGB and PNe phase, a faster wind interacts with this torus. The fast wind 
is blocked by the torus in the equatorial plane of the AGB star, but may easily flow into the 
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polar regions, creating a bipolar (P)PN. 

The development of axisymmetry, and in particular the origin of the equatorial torus 
needed to create it in the GISW model, has been attributed to various mechanisms. Here, the 
presence of a companion star, fast rotating AGB stars, magnetic fields, and the compression 
of gas into the equatorial region by two jets have been considered in various ways. We will 
briefly discuss these mechanisms below. 

A companion star may have several effects. First, it may cause the AGB star to spin- 
up (e.g. Harpaz and Soker 1994; Soker and Harpaz 2000). This spin-up may lead to small 
deviations from spherical symmetry of the AGB star itself, which are subsequently amplified 
in the dust condensation region due to the non-linear behavior and the strong temperature 
and density dependence of the dust formation process (Dorfi and Hofner 1996). This then 
leads to preferential mass loss along the equator. Second, to spin-up an AGB star, the 
companion must be close (see Sect. 2.3 for a quantification of the necessary proximity). If 
the companion is close however, and the AGB star's mass loss rate is high enough, the 
companion may accrete mass and create jets in the process. This also leads to axisymmetry 
(Soker 2005, and references therein). Finally, the gravitational influence of a companion star 
may help shape the AGB star's envelope into an axisymmetric geometry (Mastrodemos and 
Morris 1998, 1999). 

Several AGB stars and PNe show magnetic fields in their circumstellar environments, 
and these fields are often suggested to be the main agent responsible for creating axisymme- 
try. However, Soker (2006) argues that a single star can not supply the energy and angular 
momentum to create the large coherent magnetic fields required for shaping the circumstellar 
wind, although magnetic fields may have a secondary role. 

Equatorially enhanced densities are commonly attributed to the equatorially-enhanced 
mass loss of an AGB star. However, they may also originate from the compression of gas into 
the equatorial region by bipolar jets. Soker and Rappaport (2000) argue that the interaction 
of a slow AGB wind with a collimated fast wind (CFW) blown by a main-sequence or white 
dwarf companion leads to equatorial density enhancements. Here, the CFW originates from 
the accretion of the AGB wind into a disk around the companion. The CFW forms two jets 
(lobes) along the symmetry axis which compress the slow AGB wind near the equatorial 
plane, leading to the formation of a dense slowly expanding ring. Later, after the CFW and 
slow AGB wind cease, the primary star leaves the AGB and blows a second, more spherical, 
fast wind. This wind is then collimated by the dense equatorial material and leads to a 
bipolar PN as described by the GISW model. 

While these various mechanisms for axisymmetric circumstellar structures may explain 
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the shapes of PNe, they do not address the time evolution of the shape of the circumstellar 
shell. In particular, they generally fail to address the fact that most AGB shells are spherical 
(and remain that way for most of the AGB phase), whereas most post- AGB objects (and 
PNe) are bipolar (or multi-polar). 

In this paper we address the temporal changes in the axisymmetry of the circumstellar 
envelope, and in particular the development of the torus that is required in the GISW 
model. For this purpose we consider the scenario of an AGB star that rotates at a sufficient 
angular speed at the start of the AGB phase. Next, as the star evolves, we assume that the 
rotational angular momentum of the AGB star is somehow conserved. Although simplistic, 
we will demonstrate that if the above conditions are met, some very important observational 
features of AGB star axisymmtery evolution may be reproduced. In Sect. 2 we discuss our 
model and the validity of its underlying assumptions. In Sect. 3 we discuss the results of the 
model. Sect. 4 lists our conclusions. 



2. Model description 

We note that in order for what follows to succeed, the AGB star needs to have a 
sufficiently large angular velocity at the start of the AGB phase. This probably necessitates 
some sort of a spin-up, which will be discussed in Sect. 3.2. 

The mass-loss rate is affected by the escape velocity. A rotating sphere will tend to be- 
come oblate and will have a reduced gravitational acceleration, g, at the equator compared to 
the pole. With a lower g, the escape velocity is reduced and thus mass loss increases. There- 
fore we need to consider the effect this will have on the structure/shape of the circumstellar 
shell and how this effect changes with time. 



2.1. Effective gravity ratio between the equator and poles 

We start by considering a, which measures the ratio in effective gravity between the 
equator and poles of a rotating AGB star. For a test-mass m at the surface of the star, it is 
given by 

_ GMm/R 2 -mu 2 R _ u 2 R 3 
a ~ GMm/R 2 ~ l ~~GM ^' 

where G is the gravitational constant, and M and R are the mass and radius of the star. u> 
represents the angular velocity of the star at the equator. We next assume that the angular 
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momentum of the AGB star is conserved, i.e. 

Iuj = hcoi (2) 

where / represents the moment of inertia of the star. Here, we adopt the notation that x\ 
indicates the initial value of a given quantity x, i.e. the value of x right after the spin-up. 
Assuming solid-body rotation through the stellar envelope, and that at all times the star's 
deviation from spherical symmetry is small, and ignoring the core's moment of inertia and 
increasing mass during AGB evolution, we have IccM env R 2 , where M env is the envelope's 
mass. From this it then follows that 

2 2 ( -^env,i \ I R\ 



M cnv J V R 



Using Eq. 3 and 




Ti \ 2 
eff,i 



R = RAl : r\-f^) , (4) 

where L and T eff are the luminosity and effective temperature of the star, we then find for a 



a = 1 — 



i i 



GM 

where M r is the core mass of the star. 




(5) 



Using Eq. 5, we may follow a as a function of stellar evolution. Since the AGB star's 
deviations from spherical symmetry are amplified in the dust condensation region (Dorfi and 
Hofner 1996), the behavior of the degree of axisymmetry in the envelope may in principle 
be followed as a function of stellar evolution as well. 



2.2. Axisymmetry of the envelope 

The ratio p e / p p , where p e and p p represent the mass density distribution in the equatorial 
and polar direction respectively, is a measure for the degree of axisymmetry in the envelope. 
If the exact relation between p c j p v and a is known, and if we assume that p c / p p in first 
order only depends on a (see below), we may monitor p c /p P with stellar evolution. 

While the relationship between p e /p P and a is non-trivial, we can still make an estimate 
for such a relation. Dorfi and Hofner (1996) modeled the structure of a stationary dust 
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driven wind for a 1 M AGB star (T cff = 2600 K, L — 1.0 x 1O 4 L and 1.2 x 1O 4 L ) spun- 
up by a companion star to various angular velocities, up to 2 x 10~ 8 s _1 (10 _8 s~ 1 ) for the 
1.0 x 1O 4 L (1.2 x 1O 4 L ) star. Here, the spin-up leads to small deviations from spherical 
symmetry of the AGB star itself, which are subsequently amplified in the dust condensation 
region due to the non-linear behavior and the strong temperature and density dependence 
of the dust formation process. The result is an axisymmetric envelope. The models of Dorfi 
and Hofner (1996) provide rotationally-modulated solutions for the structure of stationary 
dust-driven winds, both as a function of rotation rate and polar angle. The output of their 
models includes M c /M p and v c /v p , the mass-loss rate ratio and terminal outflow- velocity 
ratio between the equatorial and polar direction. 

Using Eq. 1 and p e /p p = (M e /M v )/(v e /v v ) we calculated a and p c /p p for these models, 
in order to estimate a relation between them. The result is shown in Fig. 1. As expected, it 
can be seen that with decreasing a, i.e. stronger deviations from spherical symmetry of the 
AGB star, p c j p p increases, i.e. the envelope becomes more axisymmetric. In the remainder 
of this paper we will assume that this relation is valid for all AGB stars. The validity of this 
assumption will be discussed below. 



2.3. Validity of assumptions 

We now briefly discuss the validity of the assumption of conservation of angular momen- 
tum of the star, and the assumption of the use of the relation between a and p c j p p found in 
Fig. 1 for all AGB stars. 

The most important assumption in our model is that after the spin-up, the rotational 
angular momentum of the AGB star is conserved. The validity of this assumption is uncer- 
tain, since there are various possible sinks and sources of angular momentum for the star. 
Angular momentum will be lost from the star through its stellar wind (e.g. Soker and Harpaz 
2000). This loss can be substantial, and may in principle remove all of the stars angular 
momentum. Ignoring the magnetic influence beyond the stellar surface and the moment of 
inertia of the core relative to that of the envelop, and assuming solid body rotation through 
the envelope, Soker and Harpaz (2000) showed that the amount of angular momentum in the 
envelope is proportional to J env oc M cnv 6 , where 5 is a constant. For the upper AGB 5 = 3. 
Using the initial-final mass relation (Weidemann 2000), which links the main sequence mass 
of a low or intermediate mass star with its final mass after the AGB phase (i.e. its WD 
mass), and assuming that the AGB phase terminates when the stellar envelope's mass has 
been reduced to about 0.01 M (e.g. Blocker 1995), it can be shown that the envelope mass 
of a 1 M (7M ) changes by a factor M cnvi /M env = 45 (598), i.e. the angular momentum of 
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Fig. 1. — M c /M p (top), v e /vp (middle), and p c / p p (bottom) as a function of a. The curves 
were calculated from detailed model results by Dorfi and Hofner (1996). The calculations 
were done for a 1 M , T eS = 2600 K, AGB star for two different luminosities (for a 1.0 x 10 4 L Q 
star, rotating up to 2 x 10 _8 s _1 , and a 1.2 x 10 4 L Q star, rotating up to 10~ 8 s _1 ). For details 
see the main text. 

the star decreases by a factor J cn v,i/</cnv~ 9 x 10 4 (2 x 10 8 ) ! The loss of angular momentum 
through the stellar wind is therefore important. Other means by which the envelope may 
lose angular momentum include magnetic breaking and expansion of the convective envelope 
on the AGB (Heger and Langer 1998). 

Sources of angular momentum for the envelope may be provided by companion stars/planets 
and the core of the star. First, at the expense of its orbital angular momentum, a low mass 
companion may feed angular momentum to the envelope during a common envelope phase 
(e.g. Harpaz and Soker 1994; Soker and Harpaz 2000). The angular momentum of a com- 
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panion of mass M p is given by (Soker and Harpaz 2000) 




i 



(6) 



where Mj is Jupiter's mass and a is the initial orbital separation. Following Soker and 
Harpaz (2000), we may express the envelope's angular momentum as 



where 7 is a constant (7 = 2/9 is appropriate for the upper AGB). Comparing Eq. 6 and 
Eq. 7 shows that e.g. a M p 100 Mj = 0.1 M© companion star with an inital orbital 
separation of 1 AU may replace all of the rotational angular momentum lost by the envelope 
of a star with M ; = 1M , M cnv <O.45M , R = 1 AU, and 00 = 10~ 8 s _1 . Likewise, a 
M p 500 Mj = 0.5 M companion star with an inital orbital separation of 1 AU may replace 
all of the rotational angular momentum lost by the envelope of a star with M ; = 7M , 
M env <5.98 M , n — 1 AU, and uo = 10 _8 s _1 . The above illustrates that (a substantial 
fraction of) the angular momentum lost by the wind can in principle be replaced by a low 
mass companion star. 

Second, Garcia-Segura et al. (1999) argue that although the envelope of low- and 
intermediate-mass stars may be devoid of angular momentum (i.e. non-rotating) at the 
beginning of the thermally pulsing AGB phase, stars with main sequence masses ~1.3M 
can spin up their envelopes to rotational speeds of ~lkm/s just prior to the PN ejection. 
Here, the spin-up is caused by angular momentum from the core that leaks to the envelope 
during thermal pulses ( Garcia- Segura et al. 1999). 

The validity of our assumption of angular momentum conservation for the star thus re- 
mains to be seen. However, as we will demonstrate in the following section, if this condition is 
somehow met, some very important observational features of AGB stars may be reproduced. 
It is therefore interesting to consider the scenario where the star's angular momentum is 
conserved, despite the uncertainty of this assumption. 

Strictly speaking, the relationship found between p c j p p and a (see Fig. 1, lower panel) is 
of course only valid for the 1 M AGB star modeled above. Still, it may not be unreasonable 
to assume that p c /p p will in first order primarily depend on a, since it is the deviation from 



env 



= Ilu 



jM env R 2 uj 
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axisymmetry of the AGB star that leads to an axisymmetric envelope in the first place. If 
this assumption is correct, we may apply the obtained relationship to other stars as well, and 
generally study p c /p P as a function of stellar parameters and evolution (see below). If the 
assumption is invalid, the behavior of a still provides information on the general behavior of 
p e /pp for a given star, since deviations from spherical symmetry of the AGB star will still be 
amplified in the dust condensation region (Dorfi and Hofner 1996). Only the exact relation 
between a and p c / p p may not be specified in this case. 

3. Results of the model 

3.1. Axisymmetry and stellar parameters 

As AGB stars evolve, their luminosity increases, while their mass and effective temper- 
ature decrease. Following Eq. 5, the effect of L and T eff is thus to increase a (i.e. lower the 
degree of axisymmetry) over time, while the effect of M will be to decrease a (i.e. increase 
the degree of axisymmetry) over time. Here, the stellar mass has the largest effect on a and 
hence the axisymmetry in the envelope. While the luminosity of the star typically increases 
by a factor of 10 during AGB star evolution, a (see Eq. 5) depends on y^L\jL, which only 
decreases by a factor of 3. Likewise, while a depends on (T e s/T e ff t i) 2 , the effective tempera- 
ture typically decreases from 3500 K to 2500 K, so (T e s/T eS ^)~ 2 decreases at most by a factor 
2. In constrast, a varies with the stellar mass as (l/M)((M i — M c )/(M — M c )). Assuming its 
tip of the AGB envelope mass is ~ 0.01 M Q (Blocker 1995), a may therefore vary by a factor 
of about 3.6 x 10 3 (2.5 x 10 6 ) over the lifetime of a 1 M (7M ) star. The axisymmetry in 
the envelope will thus indeed be primarily affected by the stellar mass. 

3.2. Axisymmetry and stellar evolution 

We now investigate the time evolution of p c / p p under the influence of mass loss during 
the AGB phase. This is shown in Fig. 2 for four stars with different main sequence masses. 
It is assumed that each star is spun-up to an initial rotational velocity of 0.076 km/s. This 
choice of initial rotational velocity ensures that, when combined with the other stellar pa- 
rameters adopted for each star, a is always in the range covered by the models of Dorfi and 
Hofner (1996), therefore allowing a translation to p c /p p to be made. After the spin-up each 
star has an initial luminosity of 3000 L , an initial effective temperature of 3500 K, an initial 
radius of ~ 0.7 AU, and an initial mass loss rate of M; = 10~ 7 M Q /yr. After spin-up, the 
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Fig. 2. — p c j p p as a function of time for stars with initial masses of 1.5 M Q (dashed), 3.5 M Q 
(dashed dotted), 5.0 M (dotted), and 7.0 M (solid). For details see the main text. 

stars lose mass according to 

M(*) = M i xexp( (-L-Y\n( $ I I (8) 



where Mf is the final mass loss rate, and the time t is expressed in terms of the total AGB 
lifetime, £agb- The parameter n may be found by specifying the onset of the SW phase, tsw, 
which is also expressed in terms of £agb- Here we use 



ln(ln(^V_ln(f)) 

'sw 



- - Mj 11 V M;<- /f ,. 



where tsw is estimated from detailed model results by Vassiliadis and Wood (1993), and ^agb 
and Mf are adjusted such that the initial-final mass relation (Weidemann 2000) is satisfied 
for each star. Then the mass of the star is found by 

M(t) =Mi- f M(t')dt' . (10) 
Jo 

The exact mass-loss rate as a function of time during the AGB phase is unknown and difficult 
to determine. Still, a key feature of it seems to be that most of the mass is lost near the end 
of the AGB phase (e.g. Renzini 1981; Vassiliadis and Wood 1993). With our choice of M(t) 
this condition is satisfied. We may now insert Eq. 10 into Eq. 5 to follow a as a function of 
time. Next, we use the results of Fig. 1 to translate a into p c /p p , yielding the results shown 
in Fig. 2. Note that we kept L and T cS constant at L — lxlO 4 L and T cS = 2500 K in the 
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Table 1: Model parameters for the stars shown in Fig. 2. A£ax,2 and A£ax,5 represent the 
time before the end of the AGB phase in which p c /p P exceeds 2 and 5 respectively, i.e. when 
a reasonable axisymmetry may be observed in our model calculations. The symbol — is used 
when no value for A£ax,2 or At A x,5 is available. For an explanation of the other parameters 
see the main text. 



Mi 




L 




^eff 


Ri 




Mi 




M f 


^AGB 




A£ax,2 


At A X,5 


(M ) 


(Lo) 


(Lo) 


(K) 


(K) 


(AU) 


(km/s) 


(M /yr) 


(M /yr) 


(M /yr) 


(y r ) 


(*agb) 


(y r ) 


(y r ) 


1.5 


3000 


10000 


3500 


2500 


0.7 


0.076 


l.OxlO" 7 


l.OxlO" 5 


3.0xl0" 5 


2.3xl0 6 


0.99 






3.5 


3000 


10000 


3500 


2500 


0.7 


0.076 


1.0xl0~ 7 


1.0xl0~ 5 


3.0xl0~ 5 


2.1xl0 6 


0.96 


82 




5.0 


3000 


10000 


3500 


2500 


0.7 


0.076 


1.0xl0~ 7 


1.0xl0~ 5 


5.5xl0~ 5 


1.4xl0 6 


0.92 


202 


23 


7.0 


3000 


10000 


3500 


2500 


0.7 


0.076 


1.0xl0~ 7 


1.0xl0~ 5 


l.OxlO- 4 


6.9xl0 5 


0.80 


162 


30 



t Values estimated from detailed models results by Vassiliadis and Wood (1993). 

above analysis (where L and T ef[ are not to be confused with L ; and T cff)i ). As discussed 
earlier, the effect of L and T eS are minor compared to the effects of the stellar mass. These 
adopted values for L and T e g ensure that the derived values for p c / p p are in fact lower limits. 
Table 1 lists the parameters adopted for each model in Fig. 2. 

Fig. 2 shows that the axisymmetry increases over time. Here, the strongest increase is 
clearly only near the very end of the AGB phase, which is due to the fact that the star loses 
most of its mass only at this stage (Vassiliadis and Wood 1993; Villaver et al. 2002a). This 
result may naturally explain why most AGB stars have spherically symmetric envelopes, 
while in the post-AGB/PNe phase axisymmetry seems common-place (Meixner et al. 1999; 
Waelkens and Waters 2004; Sahai 2004; Soker and Subag 2005). 

Fig. 2 also shows that the final degree of axisymmetry, i.e. the degree of axisymmetry 
at the end of the AGB, increases with increasing main sequence mass. This result also 
seems to agree with observations. Indeed, an optical imaging survey of post-AGB stars by 
Ueta et al. (2000) with the Hubble Space Telescope (HST) seems to suggest that higher mass 
progenitor AGB stars are more likely to show bipolar structure with a completely or partially 
obscured central star than low mass stars. This result seems to be reinforced by detailed 
radiative transfer calculations of Meixner et al. (2002) for the post-AGB stars HD 161796 
and IRAS 17150-3224. 

It is interesting to note that the onset of axisymmetry does not coincide with the onset 
of the SW phase. This is illustrated in Fig. 2 and Table 1. The onset of axisymmetry only 
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becomes reasonably appearant (i.e. p c /p p ~2) in the last few tens or hundreds of years, 
long after the start of the SW phase. This is in agreement with e.g. recent findings by 
Ueta (2006). Based on Spitzer Space Telescope (SST; Werner et al. 2004) observations of 
NGC650 (Program ID 77), Ueta (2006) shows that mass loss is nearly isotropic when the 
SW is turned on towards the end of the AGB phase, and the SW mass loss precipitously 
decreases along the polar direction during the SW phase. 

The degree of axisymmetry that develops in the envelope will strongly depend on the 
initial rotational velocity of the AGB star. Generally, the rotation of single stars is expected 
to slow down significantly while they are on the Segura et al. 

1999; Soker 2006), and once these stars reach the AGB phase they are believed to rotate too 
slowly for axisymmetries to develop (Dorfi and Hofner 1996). Here, slow rotation is expected 
both with and without the presence of magnetic fields (Soker 2006). It is therefore often 
argued that a companion star is prerequisite to spin-up an AGB star (e.g. Harpaz and Soker 
1994; Soker and Harpaz 2000). The mechanism adopted in this paper would benefit from the 
presence of such a binary companion, since it may spin-up the AGB star, and be a source 
of rotational angular momentum for the AGB star's envelope that is otherwise lost through 
the AGB wind (see Sect. 2.3). Moreover, the requirement of a companion star would fit in 
nicely with the idea that binary interaction indeed plays an important role in shaping the 
axisymmetric structure of PNe (see Sect. 1), and recent observational results that strongly 
suggest that most or even all PNe are in close binary systems (De Marco et al. 2005). 

Alternatively, in case of stars with main sequence masses ~1.3M , the spin-up may 
be provided by angular momentum from the core of the AGB star during thermal pulses 
(see Sect. 2.3), in which case a companion star is not required. The thermal pulses will not 
provide the star with a spin-up at the start of the AGB phase, since these pulses only occur 
late in the AGB phase. However, they may efficiently spin-up the AGB star near the end 
of the AGB phase (Garcia-Segura et al. 1999), during which the mechanism adopted in our 
paper may still apply. In this case, it seems more appropraite to use M eilVii R3 0.1 M , since at 
the onset of the thermal pulses the mass of the envelope may already have been substantially 
reduced down to ~ 0.1 M (Garcia- Segura et al. 1999). Here, the lower initial envelope mass 
will result in a lower degree of axisymmetry in the envelope, since a (see Eq. 5) will stay 
closer to unity. 

Although the trends described for p c / p p are in nice agreement with observations, the 
absolute values predicted for p c / p p are somewhat low. For example, we predict a final value of 
Pc/Pp ~ 1-2 for a 1.5 M main sequence star. In contrast, for the post-AGB star HD 161796, 
which is believed to be of a low main sequence mass (e.g. Meixner et al. 2002, and references 
therein), Meixner et al. (2002) find p c /p P ~ 9. Also, we predict a final value of p c /p P ~ 7.3 
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for a 7.0 M main sequence star, which is on the high mass end of AGB stars. Still, Meixner 
et al. (2002) find for e.g. the high mass post- AGB star IRAS 17150-3224 p c /p p w 160, which 
is clearly much larger than predicted by our models. However, the trend of an increasing 
value of p c / p p with increasing main sequence mass is the same. 

There may be several reasons for our under-prediction of p c / p p . First, the assumed initial 
rotational velocity of 0.076 km/s is rather low. As discussed earlier, this choice ensured that, 
when combined with the other stellar parameters adopted for our stars, a is always in the 
range covered by the models of Dorfi and Hofner (1996), therefore allowing a translation to 
Pel Pp to be made. This is a computational restriction however, and not a physical one. In 
fact, rotational velocities higher than we adopted may be expected (see e.g. Garcia-Segura et 
al. 1999, and references therein). Adopting a higher initial rotational velocity will increase 
p c /pp. Also, we have only considered the effects of mass thusfar. Taking properly into 
account the effects of the luminosity and effective temperature of the star as well may also 
help increase p e /p P - For example, we adopted a luminosity of L = lxlO 4 L for all stars. 
Still, for stars with low main-sequence masses this luminosity may be too high. Lowering 
it will help increase p c / p p in this case. Finally, the values of p c / p p derived by Meixner et 
al. (2002) for HD 161796 and IRAS 17150-3224 will depend on the details of their radiative 
transfer calculations. Possibly, a different assumption of e.g. dust composition or grain 
size distribution may lower the values of p c /p P derived by Meixner et al. (2002), in better 
agreement with our results. 



4. Conclusions 

In this paper we addressed the temporal changes in the axisymmetry of AGB star 
envelopes, and in particular the development of the torus that is required in the GISW 
model. For this purpose we considered the scenario of an AGB star that rotates at a sufficient 
angular speed at the start of the AGB phase. Next, as the star evolves, we assumed that the 
rotational angular momentum of the AGB star is somehow conserved. If the above condition 
is met, some very important observational features of AGB star axisymmtery evolution may 
be reproduced. In particular compared to the star's increasing luminosity and decreasing 
effective temperature, the decreasing mass of the star primarily affects the axisymmetry of 
the envelope. Adopting a representative mass loss history, where most of the mass is lost near 
the end of the AGB phase, results in an increase of the envelope's axisymmetry over time. 
The strongest increase is only near the very end of the AGB phase, explaining why most 
AGB stars have spherically symmetric envelopes, while those of post-AGB stars and PNe 
are axisymmetric. Also, the degree of axisymmetry at the end of the AGB phase is found 
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to increase with increasing main sequence mass, and the onset of axisymmetry only occurs 
after the onset of the SW phase. Our findings are in good agreement with observations. 
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